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We propose a simple algebraic method for generating classes of traveling wave solutions for a 
variety of partial differential equations of current interest in nonlinear science. This procedure 
applies equally well to equations which may or may not be integrable. We illustrate the method 
with two distinct classes of models, one with solutions including compactons in a class of models 
inspired by the Rosenau-Hyman, Rosenau-Pikovsky and Rosenau-Hyman-Staley equations, and the 
other with solutions including peakons in a system which generalizes the Camassa-Holm, Degasperis- 
Procesi and Dullin-Gotwald-Holm equations. In both cases, we obtain new classes of solutions not 
studied before. 



PACS numbers: 02.30.Jr, 05.45.Yv, ll.10.Lm 



Finding traveling wave solutions of nonlinear partial 
differential equations has been of great interestprimar- 
ily within the context of integrable systems [l|-|4[ • Such 
studies have led to many interesting types of solutions 
in the past such as the soliton solutions, the cnoidal so- 
lutions, the compacton solutions, the peakon solutions. 
However, finding these solutions has not been easy at all 
as is evidenced in the literature. 

In a recent paper [j| , we proposed a simple method for 
generating traveling wave solutions of general nonlinear 
equations, which may or may not be integrable, starting 
from solutions of simple equations (including even linear 
equations). As we have demonstrated there through non- 
trivial examples (of mostly third order equations), this 
method is very powerful in obtaining traveling wave so- 
lutions of nonlinear equations. However, we understand 
now that the underlying reason behind the simplicity of 
our earlier proposal [f| is some very remarkable proper- 
ties of traveling wave solutions which allows us to propose 
even a simpler and more general method for generating 
classes of traveling wave solutions of a large class of par- 
tial differential equations starting from a trial traveling 
wave and an invertible map, which is the primary result 
of this letter. The earlier proposal @ forms only a part 
of this general result. 

To explain the method, let us consider a (one) space 
and time dependent nonlinear iVth order equation of the 
form 



F(v,v t ,v x ,v 2x , 



0, 



(1) 



where v(x, t) denotes the dynamical variable and the sub- 
scripts denote partial derivatives with respect to the cor- 
responding variables 



dv 



dv 



d n v 



Vnt- 



<9™ H 



dt n dx m 



(2) 



We note that equation {1} may or may not be related 
to an integrable system. A traveling wave solution of 
the system has the form v(x, t) = v(kx — uit), where k 



denotes the wave number and uj the frequency of the 
traveling wave with the dispersion relation determined 
by the dynamical equation. 

We note that for any given form of the traveling wave 
solution v{kx — uit), all the x and t derivatives can be 
expressed in terms of the solution itself. Thus, for a 
traveling wave we can write 



v x = Vl(v), vt = -{u/k)Vi(v), 



(3) 



where the functional form of V\ (v) can be explicitly deter- 
mined if the form of v is known. Let us denote y = kx—ut 
and explain this assertion a bit more in detail. If v = f(y) 
is defined by a (differentiable) bijective map (one-to-one 
and onto) then / has an inverse and y — In this 



case, we can write 



dv 



^=vy = nr\v))= g (v). 



(4) 



However, as we know many maps (/) in physics are not 
invertible @. If / is not a bijection (so that its inverse 
does not exist), then we can still define an inverse of / 
on each domain Di where the function /(y) is monotonic 
(so that an inverse f~ x exists on the domain Di). In this 
case we can write 

^=/'(/r 1 («))=5i(«), for 2/ € A, (5) 

where fT~ is the inverse of / on the domain Di. The 
assertion in (j3|) makes this assumption. Furthermore, 
if we make the identification v = Vq, then any higher 
derivative of v can be written as 



V-at,; 



v, 



(G) 



where for continuous functions, V n (v) are related to V\(v) 
recursively as 



V n (v) = Vi(v), 

dv 



n > 1. 



(7) 
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For example, we note that 

dv x dVi(v) dVi(v) dV t (v) ... . 

V 2 = v xx = — = — = — v x = — Vi(v), 

ax ax av av 

. ® 

and as a result, the expression ([7J can be used recursively 
to express V n as a function of v when the form of V\ (y) 
is known. 

Therefore, we see that for a traveling wave solution 
v(kx-Lot), the dynamical equation (fl} of order N reduces 
to an algebraic equation of the form 



F(u),k;v,V u V 2 ,-~ ,V N ) = 0. 



(9) 



This equation leads to relations between w, k and the 
functions V n (v) and can also determine the dispersion re- 
lation for a traveling wave solution. The algebraic equa- 
tion (jSJ) can be the starting point to check if a given 
ansatz v(kx — cot) solves ([1]). Namely, with a given ansatz 
for v, the form of V\{v) and, therefore, all the V n (v) would 
be determined and ([9]) would reduce to a polynomial al- 
gebraic equation in v which is undoubtedly much easier 
to analyze. However, our goal is to generate classes of 
traveling wave solutions with parameters which for dif- 
ferent values lead to different kinds of solutions and we 
propose to do this in the following manner. 

The method basically involves two steps. First, we 
choose a simple traveling wave u(kx — Lot) which we call 
a trial wave. Then we introduce an invertible map 0] 



u(kx — Lot) = J-{v(kx — Lot)), 



(10) 



relating the trial wave to the kind of traveling wave solu- 
tion we are looking for. It is worth noting here that this 
procedure may seem rather arbitrary. However, once a 
trial function is chosen (trigonometric, hyperbolic, • • • ) 
which is simple to manipulate and we keep in mind the 
kind of solution (soliton, compacton, peakon, • • • ) of the 
dynamical equation that we are interested in, the map 
can be constructed systematically with little arbitrari- 
ness (as we will describe in the examples below). The 
map is really chosen keeping the nature of the solution 
of interest in mind and may involve several parameters. 
Requiring v(kx — Lot) to satisfy the dynamical equation 
([1]) or equivalently ([9|) may determine some of these pa- 
rameters while the undetermined parameters would lead 
to a class of solutions depending on these parameters. 

Checking if v{kx — Lot) given by the inverse of the map 
in (llOj) satisfies the dynamical equation is quite easy in 
its algebraic form ([§]). For example, the trial wave will 
lead to a relation of the form u x = U\(u) (see ([3])), where 
Ui(u) is known explicitly from the known form of the 
trial wave and the map (fTO)) leads to 



r{v) ■ 



d7» 
dv 



(11) 



Namely, the known from of the trial wave u as well as the 
choice of the map T would determine the (polynomial) 
forms of all the unknown functions V n (v) upon using ([7)1, 



As a result, the dynamical equation (J9j) will reduce to 
an algebraic (polynomial) equation in v. The solution 
of this simple algebraic equation would possibly deter- 
mine some (or all) of the parameters in the map and 
would generate a traveling wave solution of the original 
dynamical system ([I} through the inverse map (see ([TO]) ) 
v(kx — Lot) = J 7 [u{kx — Lot)). If all the parameters of 
the map are determined by the equation, then we have 
a particular traveling solution of the system (as in [5j). 
However, if the map is chosen to be rich enough and some 
of the parameters are left undetermined by the equation, 
then this can lead to interesting classes of traveling wave 
solutions as we will illustrate through examples. 

It is clear that this method applies equally well to equa- 
tions which may or may not be integrable, and so the 
potential for its applicability is indeed quite great. How- 
ever, we point out here that if there are discontinuities 
in the solutions (as some of the interesting solutions do) , 
then the simple recursion relation J7J may not hold and 
a direct evaluation of the higher derivatives may become 
necessary. This would be true in general when the deriva- 
tive v x (or u x ), which leads to V± (or f/i), is not a single- 
valued function of v (or u). On the other hand, if the dy- 
namical equation (JlJ does not involve v x directly, rather 
f(v x ) which is free from ambiguity, then our method can 
be carried out easily as we explicitly show in the examples 
below. 

Let us explain the method as well as its potential with 
a couple of examples. We first study traveling waves 
of the form of compactons which are configurations of 
finite extent free of exponential tails and which arise as 
solutions of the Rosenau-Hyman (RH) equation Q . Here 
we would study the deformed equation 



v t + (P{v)) x + {Q{v)v xx ) x = 



(12) 



where P(v) — a\V + a 2 v 2 + o 3 v 2+p + v 2+2p , and Q(v) — 
v l+2p , with ai, 0:2,0:3, and p > real constant param- 
eters. Equation (p~2|) is motivated by the RH family of 
equations v t + (v m ) x + (v n ) xxx = Q, as well as by 
the N dimensional Rosenau-Pikovski (RP) @ and the 
Rosenau-Hyman-Staley (RHS) [1(3] family of equations 
v t + {v m ) x + (l/s)(v r V N v s ) x = 0, where a unidirec- 
tional convection is balanced by a ./V-dimcnsional dis- 
persive force. Although the deformed model (|12p is one 
dimensional, we will show below that the nonlinear con- 
vection and dispersion which are present in (1121) bring 
a variety of effects as we change the parameters. Since 
the results are analytical, our procedure may be of some 
use in the search of analytic solutions in higher dimen- 
sions, an issue which is beyond the scope of the present 
work. For a traveling wave solution of (fT2"|) . we proceed 
as discussed earlier and obtain the relation (see ©) 

I - P'{v) - Q'(v)V 2 (v) - Q{v)V 2 \v) = 0, (13) 

where a prime denotes derivative with respect to the ar- 
gument as in dill) . 
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To study compactons, let us choose the trial wave in 
the form u{kx— wt) = cos(kx—wt), for \kx— wt\ < it, and 
— 1 otherwise, which leads to u x — U\(u) = — fcsgn(fcx — 
wt) \/l — u 2 for \kx — wt\ < 7r. The presence of the sgn 
function seems to introduce an ambiguity, but we note 
that (fT5|) only depends on V 2 and its derivative so that 
there is no ambiguity in the equation. 

Let us next introduce the map u = a + v p , where a,p 
are real constant parameters so that we can write v{kx — 
Lot) = (cos{kx — ust) — a) 1 / p . In this case, for \kx — wt\ < tt 
the trial wave and its derivatives are continuous and we 
can follow our general procedure. Thus, we use ([7]) and 
(fTTj) to determine 

V 2 (v) = ^v + a{p-2)v 1 -P+(l-p)(l-a 2 )v 1 - 2 P). (14) 
Using this in (|13[) we obtain the algebraic equation 
(a 1 -^)+2^ 2 -^(l-p)(a 2 -l)^ 
+ (2 + p) (a 3 -^(2-p)ay i+ P 

+2(1 +p) (l-^J)« 1+2p = 0. (15) 

This determines the parameters as k = ±p, w = ±aip, 
a 2 = (1 — p)(a 2 — 1), and a 3 — (2 — p)a. 

We can now construct a variety of solutions with dif- 
ferent choices of the above parameters. Let us also note 
here that since a constant is a solution of (fTS"]) . by adjoin- 
ing such a solution to the solution found with the inverse 
map in a continuous manner, one can easily construct 
new solutions of the dynamical system. For example, 
with p an odd integer or the inverse of an odd integer, 
we obtain the compacton solution 

{(cosp(x~ ait) ~ a)p, for \p(x — ait)\ < tt, 
— (1 + a) p , otherwise, 

which satisfies (fl"2j) . Note that if we redefine v — >• v + b 
in (fl"2"|) . this new parameter b can be used to control the 
asymptotic behavior of the solution, such that for b = 
(1 + a) 1 /? the shifted solution will vanish outside the 
interval [— tt, tt]. 




FIG. 1: Plots of the traveling wave solution for p = 1/3, with 
a — 0, a — 0.3 and a = 0.8 (left panel), and with a = —0.2, 
a = —0.4 and a = —0.8 (right panel). 

There is a variety of traveling wave solutions contained 
in (|16[) and some of them lead to interesting features 



which we now illustrate. In Fig. 1 we have plotted the 
solution (|16p at t = for p = 1/3, and for some values of 
a G [—1,1], where the amplitude increases with increas- 
ing \a\. We note that in the limit a — > — 1 + we have a 
compacton of standard bell shape. As a increases, the 
solution develops a lump at the top which is a nice nov- 
elty for compactons, with the interesting feature that the 
height of the bottom (hi,) and top (h t ) portions of the 
new compactons obey: h\, < h t for a G (—1,0), hi, = h t 
at a = 0, and hi, > h t for a £ (0, 1). We further note that 
in the limit a — > 1~, the solution leads to a compacton 
with a flat plateau. 

In Fig. 2 we have plotted solutions at t = for p = 3, 
and for some values of a, and there the amplitude in- 
creases with decreasing \a\. These plots unveil remark- 
able features of the traveling waves: for a = ±0.5 the so- 
lution is a [thin (+) or thick (— )] tipon-like compact wave 
Q, for a — 1 it is a peakon-like compact wave [ll| (see 
below), and for a = —1 it has an oval- like compact shape. 
All of the above solutions identify distinct types of com- 
pactons, some of them have never been studied before. 
In both cases, we note that the constant a, introduced 
by the chosen map, works like a deformationparameter, 
since it deforms the shape of the solution Q . 




FIG. 2: Plots of the traveling wave solution for p — 3 and for 
a = 0.5, a = 1 and a = 1.5 (left panel), and for a = —0.5, 
a = —1 and a = —1.5 (right panel). 

To introduce the second example, let us recall that the 
equation v t - v xxt + (r + l)vv x - rv x v xx - vv xxx = is 
known to be integrable for r — 2,3. For r = 2 it is known 
as the Camassa-Holm (CH) equation |ll| which was de- 
rived as an approximation to Euler's incompressible fluid 
motion. For r = 3, it describes the Degasperis-Procesi 
(DP) equation [l2[ and both these equations represent 
models of shallow water waves Furthermore, both 
these equations are known to possess peakon solutions 
which have discontinuous first order derivatives. Since 
solitary waves constitute a well known form of the wa- 
ter waves, another integrable equation was recently in- 
troduced by Dullin-Gattwald-Holm (DGH) [l3j], which 
contains the KdV equation as well as the CH equation in 
different limits, thereby having the interesting feature of 
supporting both the bell shaped solution of KdV as well 
as the peakon solution of CH as limiting cases. 

Keeping this in mind, let us study, as another example 
of our method, a deformed system of the form 

vt - v xxt + f(v)v x - g{v)v x v xx - h{v)v xxx = 0, (17) 
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with 

f(v) =/3 +/3 1 v + /3 2 v 2 , 

g(v)=r + sv, h(v) = p + Piv + /3 2 v 2 . (18) 

We note that when /3o = 0o = P% = P2 = s = 0, /3i = 
(r + 1) and /?i = 1, equation (fl~7|) reduces to the CH-DP 
equations (for r — 2,3) while for $2 — P2 — s = it leads 
to the DGH equation [l3j]. For general = 0, 1,2 

and r, s the deformed equation (jTTJ) is indeed a much 
richer system and to clarify our method we try to find 
the traveling wave solution of this generalized system. 

A traveling wave solution of (fl~7| of the form v(kx— u>t), 
would satisfy a relation of the form (see © ) 



ui 



'0 - 0iv - /3 2 v 2 + (r + sv)V 2 
+ ( A> + Piv + f3 2 v 2 - |) V{ = 0. (19) 



Keeping a soliton solution in mind, let us next choose a 
simple trial wave of the form u(kx — Lot) — cosh(fca; — tut) . 
Here we note that only V 2 and its derivative are present 
in (|19p so that the apparent ambiguity with the sign in V± 
is again not present. With this, we choose a map (|10p of 
the form (keeping the soliton in mind) u = \fajv, where 
a is a constant parameter to be determined. In this case, 
following our earlier discussion, we obtain 



V 2 = Ak 2 v - 6k 2 v 2 /a 



(20) 



As a result, (|T9|) leads to the algebraic equation for v 
of the form 



i-*)-«(i-sO*) 

+ ^12(|~^o)^ + 4(r + ^)fc 2 -/?^i; (21) 

/ / „ „ v fc 2 \ ~ k 2 

- f/32+(6(r+2/3i)-4a(s+/3 2 )j — Jw 2 -6(s+2/3 2 )— v 3 =0. 



Requiring the coefficients of each power of v to vanish 
now determines the parameters as well as the traveling 
wave solution for ([TT]). 

There are several solutions of (|2"Tj) depending on the 

values of the parameters in (fT7)l . If /3o = /?o, fii = s = 
and j3 2 = r + /3i ^ 0, then fl^T) determines k = ±1/2 and 

w=± — — ± y a= — 2^' (22) 

and the traveling wave solution corresponds to 



2f3 2 



■sech 2 \ -qp 



4/3 2 



(23) 

We note from this that there are several distinct choices 
of parameters which lead to waves of both positive and 



negative velocities. Furthermore, in our approach, fi 2 is 
merely a deformation parameter and if we choose f3 2 < 0, 
we indeed have a bell shaped solution. 

On the other hand, for = s = 0, u = kj3o — k(3o we 
have 



±- 



1 



/ & 3(r + 2/3 1 )/3 1 

~ 2 '\r + Pi Q 2(r + j8i)j8a 

with the traveling wave solution 



(24) 



01 



•(x-fat). (25) 



= _3(r + 2W lsech2 l 
2(r + 0i)02 2 



We note that the velocity in this case is determined by 
the constant /3o alone. There are other choices of param- 
eters, leading to distinct bell shaped solutions, but we 
now focus attention on the peakon solution. 

To obtain the peakon solution for the deformed equa- 
tion (fl7|) . let us choose the trial wave to have the form 
u(kx — cot) = exp(— \kx — u)t\), which leads to u x = U± = 
—ksgn(kx-wt) u. We note that the solution has discon- 
tinuous derivatives and as mentioned earlier, the recur- 
sion relation does not hold. In this case, we need to 
calculate the higher functions explicitly. Let us rewrite 
(jTTJ) in the form 

((/(«) - |) - 9(v)V 2 ) Vr = (h(v) - |) V 3 , (26) 

and note that the presence of V\ in the left hand side 
and of V3 in the right hand of this equation removes any 
ambiguity with the sgn function. The next step is to use 
(fTTj) and u nx = U n (u) to calculate 



V 2 (v) = U 2 {F{v))/T' ~V 2 F"/F', 
V 3 (v) = U 3 (T(v))/T' -W X {V2F' + V 2 F") 2 F" 
+V 1 3 (3T" 2 -T'T"')/T' 2 . 



(27a) 
(27b) 



The higher derivatives can also be calculated easily, but 
they are not necessary in the present case. We then 
choose the map (|10p u = (v — b)/a, where a and b are 
real parameters. With this, after substituting the explicit 
forms of Vi , V2, V3 we end up with the algebraic equation 

*sga(te-wt)(«-6)n (~A))-(|-A) + r6) k 2 ^v 

+((r + ft - sb)k 2 - p x )v 2 + ((a + p 2 )k 2 - fo)v^j 

+2k 3 aS'(kx - u)t) (j3 + Av + /3 2 v 2 - ^ =0. (28) 

Here we have simplified some of the terms using the con- 
ventional relation sgn(fca; — Lot)8{kx — cot) = 0. Also, 
using the identity kS'(kx — oot)v n — (S(kx — oot)v n ) x — 
k(a + b) n 5'(kx — cot), the last term in (|28|) can be simpli- 
fied even further. 

There are distinct solutions of (|28p depending on the 
values of the parameters in (|17p , but the last term in 
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determines u> = k((3 + Pi(a + b) + I3 2 (a + b) 2 ). We note 

that for b — 0, and for /3 = /?o, /?i = r + /3 2 = s + flz, 
equation (f2"5| leads to A: = ±1 without any restriction on 
a and the traveling wave solution has the form 

v(x,t) = ae-\ x ~ U3t y (29) 

If we further take the simple choice (among many) of 

the parameters, /3o = @2 = 0,/?i = ±1, we can identify 
a = uj and the solution (I29p coincides in this case with 
the peakon solution in 11,]. However, a nice novelty of 
our solution (j2"9")l is that the amplitude does not need to 
be equal to the velocity anymore, and as a result, we 
can also have a peakon solution with positive amplitude 
traveling with positive or negative velocity. For 6^0, 
there are peakons where b can be used to control the 
width and the asymptotic behavior of the solution. All 
the above results show a rich diversity of solutions of 
(fTT)) some of which, to the best of our knowledge, have 
not been studied in the literature before. 

In conclusion, in this letter we have proposed a simple 
method for constructing classes of traveling wave solu- 



tions of nonlinear partial differential equations. Exploit- 
ing the properties of traveling waves, the method converts 
the dynamical equation into a simple algebraic equation 
which is easy to analyze. The two steps of taking a simple 
trial wave and then choosing a map (with several param- 
eters) between the trial wave and the unknown traveling 
wave solution lead in a very efficient way to classes of 
traveling wave solutions. This method is a generalization 
of our earlier results, and we have demonstrated the ver- 
satility and the potential of our method by constructing 
the compacton solution for the Rosenau-Hyman equation 
where we have also shown the existence of newer kinds of 
compacton solutions, and the bell shaped solution as well 
as the peakon solution for the Camassa-Holm equation, 
bringing novelties for both the bell shaped and peakon- 
like solutions. In particular, we have found some very 
interesting new compact solutions, which are the tipon-, 
peakon- and oval-like compactons. The method will be 
further used to study other nontrivial ex amp les of current 
interest in nonlinear science separately (l4| . 

The authors thank CNPQ, PRONEX-CNPq-FAPESQ, 
and the US DOE Grant number DE-FG 02-91ER40685. 
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